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plate  into  a  uniform  crossing  flow  was  formulated  to  provide  a 
simplified  method  of  predicting  the  interference  effects  arising 
from  the  complex  flow  fields  induced  by  ship  bow  thrusters.  This 
model  was  an  extension  of  previous  work  based  upon  a  description 
of  the  jet  as  a  series  of  distributed  vortices.  The  analysis 
takes  into  account  the  position  of  the  effective  source  of  the 
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crossflow.  For  representative  jet-to-crossflow  velocity  ratios, 
the  flow  field  and  pressure  distributions  were  calculated 
utilizing  different  combinations  of  effective  source  nosition 
and  blockage.  The  aacuracy  of  the  model  was  evaluated  by  com¬ 
parison  with  the  available  experimental  data.  Although  good 
agreement  was  achieved  for  large  portions  of  the  interaction 
field,  several  regions  were  identified  requiring  further 
analytical  description. 
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ABSTRACT 


An  analytical  model  of  a  jet  injected  normally  from  a  flat 
plate  into  a  uniform  crossing  flow  was  formulated  to  provide 
a  simplified  method  of  predicting  the  interference  effects 
arising  from  the  complex  flow  fields  induced  by  ship  bow 
thrusters.  This  model  was  an  extension  of  previous  work 
based  upon  a  description  of  the  jet  as  a  series  of  distributed 
vortices.  The  analysis  takes  into  account  the  position  of 
the  effective  source  of  the  jet  and  the  blockage  due  to  the 
presence  of  the  jet  in  the  crossflow.  For  representative 
jet-to-crossflow  velocity  ratios,  the  flow  field  and  pressure 
distributions  were  calculated  utilizing  different  combinations 
of  effective  source  position  and  blockage.  The  accuracy  of 
the  model  was  evaluated  by  comparison  with  the  available 

i 

experimental  data.  Although  good  agreement  was  achieved  for 
large  portions  of  the  interaction  field,  several  regions  were 
identified  requiring  further  analytical  description. 
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I.  INTRODUCTION 


In  recent  years  there  has  been  an  increase  in  the  instal¬ 
lation  and  use  of  hull  mounted,  tunneled  bow  thrusters  to 
improve  ship  maneuverability  in  restricted  situations,  such 
as  dynamic  position  holding  and  confined  waters  docking. 

This  increase  has  created  renewed  interest  in  the  development 
of  am  ability  to  predict  bow  thruster  performance.  To  date, 
there  has  been  relatively  little  pre-installation  design 
consideration  regarding  precise  power  requirements  or  ultimate 
operational  performance.  The  limited  testing  that  has  been 
conducted  /"“l-£7*  has  resulted  in  a  recurring  major  problem, 
namely,  radical  variations  in  effective  turning  moment  with 
minor  changes  in  ship  speed.  Chislett  and  Bjorheden  /~4  7 
noted  that  "an  area  of  low  pressure  is  created  downstream 
of  the  discharging  jet  producing  a  resultant  suction  force 
that  has  a  shifting  center  of  action  with  changes  in  ship 
speed,”  thereby  altering  the  effective  moment  produced  by 
the  thruster.  The  ability  to  predict  the  magnitude  of  this 
resultant  suction  force  and  its  center  of  action  would  permit 
designers  to  accurately  determine  powering  and  control  re¬ 
quirements  for  thrusters. 

This  problem  of  radical  variations  in  effective  side 
forces  and  turning  moments  with  only  minor  changes  in  ship 
speed  may  be  better  understood  by  considering  the  general 

1 - 

Numbers  in  brackets  refer  to  references  listed. 
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equations  of  motion  of  a  ship  moving  in  the  horizontal 

plane  of  the  sea.  From  /“5  7  these  equations  ares 

m(u  -  rv)  *  X 

m(v  -  ru)  =  Y 

I  r  *  N 
z 

where  m  -  mass  of  the  ship 

u  ■  velocity  in  the  x-direction 

v  *  velocity  in  the  y-direction 

r  »  angular  velocity 

I  *  mass  moment  of  inertia 

X  *  excitation  force (s)  in  x-direction 

Y  *  excitation  force (s)  in  y-direction 

N  =  excitation  torque (s) 

(dots  indicate  differentiation  with  respect  to  time) 

To  minimize  the  complexity  of  the  following  brief  analysis, 

l 

all  forces  and  torques  will  be  referred  to  the  ship's  center 
of  gravity.  See  Fig.  1.  < 

When  these  equations  are  transformed  into  their  linearized 
version  /"“  5_7>  they  become 

~  xu(u  "  Uo*  +  (m  "  Xu)  uo  ■  X 

-  Yyv  +  (m  -  Y^)v  -  (Yr  -  muQ)r  -  Y*r  -  Y 

-  Nvv  -  N^v  -  Nrr  +  (I2  -  N»)r  -  N 

•  •  • 

where  the  subscripts  u,  u,  v,  v,  r,  r,  denote  differentiation 
with  respect  to  the  variable  subscripts,  and  uQ  is  constant 
reference  surge  velocity.  The  left  sides  of  these  equations 
explicitly  represent  the  velocity  and  acceleration  dependent 
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reactions  due  to  drag,  added  mass,  added  inertia  and  other 
coupled  motions;  while  the  right  sides  implicitly  represent 
excitation  force (s)  and  moment (s) ,  which  are  in  this  analysis 
due  to  the  thruster.  Since  this  is  a  linearized  formulation 
utilizing  small  angle  perturbations,  the  excitation  force  in 
the  x-direction  is  negligible  compared  to  that  in  the 
y-direction: 


X  <<  Y 

therefore  X=0.  The  y-direction  excitation  force  can  be 
expressed  as 

Y  *  Yj  +  Y±  (1) 

where  *.’j  ■  itkL\,  the  ideal  jet  thrust  produced  by  the 
thruster,  and  Yi  is  the  induced  force  created  by  the  inter¬ 
action  of  the  thruster  discharge  with  the  flow  past  the  ship. 
Utilizing  experimental  data  from  /"” 4_7  in  the  range, 

0.1  i  <  0.25,  Y.  can  be  approximated  by  a  linear  function 

j 

of  ship  speed 


0.22  +  k(g-) 


(2) 


where  -3.0  <  k  <  -1.5.  Substituting  a  nominal  value  of 
k  »  -2.0  and  Yj  *  ^nto  equation  (2)  gives 

Yi  -  mjDj[o.22  -  2.0  <§-)] 


(3) 


Substituting  equation  (3)  into  equation  (1)  results  in 


Y  ■  m.U  .  [  1.  22  -  2.0(?r-)  1 


which  explicitly  couples  the  y-direction  motion  to  the 
x-direction  motion,  normally  uncoupled  in  the  linear 
approximation.  In  addition  to  the  coupling,  in  the  relevant 
velocity  ratio  range,  0.1  <  jf-  <  0.25,  there  is  an  appreciable 
reduction  in  the  net  thruster  force  due  to  motion  in  the 
x-direction.  Furthermore,  the  moment  produced  by  the  thruster 
is  a  function  of  the  net  effective  thrust,  Y,  and  a  character¬ 
istic  moment  arm,  xfc: 

Nt  -  XfcY 

As  reported  by  Ref s .  /~  1-3/,  the  degradation  of  net  effective 
thrust  is  accompanied  by  a  shift  in  the  center  of  action  of 
the  induced  force,  Y^,  which  results  in  a  change  in  the  char¬ 
acteristic  moment  arm.  This  shift  in  the  thruster's  net 
effective  thrust  center  of  action  is  a  direct  consequence  of 
the  coupling  of  x-direction  motion  with  turning  moment. 

This  simplified  analysis,  which  neglected  all  excitations 
other  than  the  thruster,  points  out  the  complex,  coupled  mo¬ 
tions  arising  from  the  interaction  of  the  thruster  discharge 
and  the  flow  past  the  ship.  To  fully  understand  this  highly 
non-linear,  three-dimensional  flow  situation  extensive 
analytical  and  experimental  work  is  required.  This  thesis 
will  attempt  to  provide  some  understanding  of  the  interaction 
of  a  jet  discharging  into  a  uniform  cross-flow  in  order  to 
develop  some  ability  to  predict  the  degradation  of  net  effective 


thrust 


II.  BACKGROUND 


A  survey  of  the  literature  pertaining  to  the  installation 
and  operation  of  bow  thrusters  revealed  a  surprising  lack  of 
experimental  data,  considering  the  number  of  units  currently 
installed.  English  /T,27#  Stuntz  and  Taylor  /”3 _7,  and 
Chislett  and  Bjorheden  /”4 7  present  some  data  from  model  tests 
conducted  at  various  facilities,  but  in  these  works  there  is 
no  formalized  recording  of  full  scale  test  or  actual  instal¬ 
lation  test  data.  However,  further  investigation  led  to  a 
related  field:  Vertical  or  Short  Take-Off  and  Landing  Aircraft 
(V/STOL)  research.  The  concept  of  using  air  jets  positioned 
in  aircraft  wings  to  provide  vertical  take-off  and  landing 
capability  correlates  extremely  well  with  flow  patterns  pro¬ 
duced  by  a  bow  thruster.  As  reported  by  Bradbury  and  Wood 
Ce  7  for  incompressible  flow,  the  jet  (thruster  discharge) 

path  and  induced  flow  are  dependent  mainly  on  the  momentum 
2  2 

flux  ratio,  PjUj  /pu  ,  and  independent  of  Reynolds  number. 
Margason  /"" 7_7  and  Gordier  /"" 8_7  also  reported  that  the 
effective  velocity  ratio  is  the  predominate  characteristic 
in  determining  the  path  of  the  jet.  Therefore,  until  more 
extensive  ship  or  model  data  have  been  collected,  the  re¬ 
sults  of  V/STOL  research  are  used  as  a  basis  for  predictive 
theories. 

Although  there  has  been  a  good  deal  of  research  conducted 
in  the  area  of  air  jets  in  cross-flows,  the  majority  of  the 
work  has  been  concerned  with  analytically  or  empirically 
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defining  the  jet  trajectory  and  jet  cross-section  geometry. 
The  interference  effects  of  the  interacting  flows  have  re¬ 
ceived  relatively  little  attention.  A  fluid  jet  injected 
into  a  crossing  stream  has  previously  been  described  as  a 
turbulent,  three-dimensional,  highly  non-linear  flow  problem 
which,  even  with  considerable  simplification,  requires  ex¬ 
tensive  computer  time  to  numerically  solve  the  appropriate 
form  of  the  Navier-Stokes  equations  /—9 _7.  Therefore,  early 
investigators  have  sought  simplified  methods  for  determining 
the  interaction  effects,  such  as  pressure  and  velocity  dis¬ 
tributions,  and  have  realized  the  necessity  of  knowing  the 
geometry  of  the  jets  trajectory  and  shape. 

Experimental  data  form  the  basis  for  most  trajectory 
formulations.  Either  pressure  or  velocity  measurements  have 
been  taken  in  the  flow  field  encompassing  the  jet-crossflow 
interaction  region  and  curves  fit  through  points  of  maximum 
pressure  or  velocity.  These  curves  were  then  compared  with 
flow  visualizations.  In  some  instances,  semi-empirical  for¬ 
mulations  were  derived  from  the  conservation  laws  combined 
with  experimentally  determined  constants.  Abramovich  /""  107 
qualitatively  described  the  turbulent  jet  in  a  deflecting 
flow  and  presented  some  empirical  methods  for  predicting 
the  trajectory  of  the  deflected  jet.  Jordinson  /Tl  7  pre¬ 
sented  trajectory  data  by  recording  and  plotting  contours 
of  total  pressure  coefficients.  Keffer  and  Baines  ^”12  7 
presented  experimental  results  from  which  it  was  determined 
that  for  various  jet  strengths,  the  jet  trajectories  could 
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be  represented  by  a  single  function.  Gordier  /~ 3_7,  the 
only  reported  source  of  water  jet-water  crossflow  experiments, 
presented  an  empirical  jet  trajectory  formulation  based  on 
a  curve  passed  through  experimental  points  of  maximum  total 
pressure.  Sucec  and  Bowley  /“9  7  formulated  an  analytical 
expression  for  the  jet  trajectory  utilizing  previous  experi¬ 
mental  information  and  the  assumption  that  the  distributed 
pressure  force  and  entrained  momentum  flux  could  be  approxi¬ 
mated  by  an  aerodynamic  drag  force.  Margason  /”7 _7  utilized 
flow  visualization  techniques  and  pressure  measurements  to 
formulate  a  trajectory  equation.  In  addition,  his  paper  also 
reviewed  and  compared  the  results  of  a  number  of  other 
studies  concerned  with  trajectory  prediction.  All  the  above 
predictive  formulations  produced  trajectories  within  a  range 
of  uncertainty  that  can  be  attributed  to  expected  experimental 
error,  differences  in  test  procedures  and  facilities,  and  in 
the  case  of  analytical  formulations,  simplifying  assumptions. 
In  view  of  this,  it  was  concluded  that  any  one  of  them  re¬ 
flected  the  current  state  of  the  art  in  trajectory  prediction. 

The  ability  to  predict  the  jet  trajectory  is  of  definite 
importance,  but  the  primary  objective  of  this  jet-crossflow 
research  is  to  provide  insight  into  the  changes,  in  the  other¬ 
wise  uniform  flow,  created  by  injecting  a  jet.  Experiments 
conducted  by  Bradbury  and  Wood  /“ 6_7,  Vogler  £“13 _7,  Fearn 
and  Weston  £“14_7,  McMahon  and  Mosher  15_7,  Kamotani  and 
Greber  ^”16  7  and  others,  provide  measured  values  of  the 
pressure  distribution  on  the  surface  surrounding  the  jet 


orifice  which  have  been  used  to  compare  the  accuracy  of 
predictive  models.  During  the  past  few  years ,  a  number  of 
approximate  methods  to  predict  the  pressure  distribution  has 
been  developed:  Wooler  /T7,18 7,  Wu  and  Wright  /T97,  Adler 
and  Baron  /”20 7  and  Schmitt  /"*: 2 17*  Again,  because  of  the  com¬ 
plex  nature  of  the  flow,  all  these  models  have  utilized  some 
degree  of  empiricism.  Either  the  model  is  formulated  using 
experimentally  obtained  trajectory  expressions  or  the  model 
is  formulated  using  analytical  expressions  whose  coefficients 
are  selected  to  provide  best  fit  with  experimental  data. 

The  majority  of  the  models  are  based  on  integral  techniques 
with  simplifying  assumptions,  such  as  (i)  the  representation 
of  the  entire  flow  as  two-dimensional  /T9 7,  (ii)  the  external 
flow  is  irrotational,  uncompressible  and  steady-state  /~20 7  or 
(iii)  the  flow  is  inviscid  except  that  viscosity  is  the 
mechanism  that  leads  to  entrainment  /T8 _7. 

With  the  aid  of  these  assumptions,  the  following  approaches 
were  made.  Adler  and  Baron  formulated  their  model  by  integrat¬ 
ing  momentum  equations  to  describe  the  jet  mixing  field,  without 
using  empirical  trajectory  data.  However,  numerous  parameters 
were  derived  using  empirical  correlations  from  prior  research. 
Their  results  provide  satisfactory  agreement  with  experiment, 
but  from  the  quantity  of  empiricism  involved,  one  might  ex¬ 
pect  such  agreement.  Wooler  £“1  lj  followed  by  Wu  and  Wright 
£~1 9_7  utilized  blockage-sink  representations  where  entrain¬ 
ment  of  crossflow  fluid  was  handled  analytically  with  suitably 
chosen  coefficients.  Again  the  results  had  satisfactory 
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agreement  with  experiment.  However,  an  earlier  effort  by 
Vfooler  /~18_7  minimized  empiricism  to  the  extent  that  only 
an  experimentally  determined  jet  trajectory  was  necessary 
to  complete  the  model.  By  thus  restricting  the  number  of 
empirical  parameters,  this  model  has  considerable  appeal 
from  an  engineering  point  of  view.  Since  visualization  of 
actual  jet-crossflow  situations  has  shown  that  the  jet  is 
deflected  in  the  crossflow  direction  and  forms  two  contra¬ 
rotating  trailing  vortices,  this  interference  model  was 
based  on  the  representation  of  the  jet  by  a  distribution  of 
vorticity.  Incorporating  arguments  from  the  aerodynamic 
theory  of  lifting  bodies  this  vorticity  distribution  was 
quantitatively  described  along  the  experimentally  determined 
trajectory.  By  virtue  of  minimal  empiricism  and  therefore 
greater  appeal,  this  model  was  selected  as  the  basis  for 
continued  investigation  in  developing  a  better  predictive 
tool. 
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III.  REVIEW  OF  WOOLER'S  VORTEX  MODEL 

Before  any  attempt  to  develop  an  improved  model  was 
made,  a  detailed  re-derivation  of  the  basis  formulation  was 
performed  in  order  to:  (1)  gain  a  better  understanding  of 
the  logic  of  formulation  and  (2)  confirm  the  reported  results. 
A  review  of  the  distributed  vortex  method  follows. 

One  area  of  interest  was  the  geometry  of  the  problem. 

The  empirical  trajectory  equation  used  was 

x/d  =*  B  £cosh(z/Bd)  -  l]  (4) 

where  x  is  measured  in  the  direction  of  the  mainstream  and 

z  is  measured  in  the  direction  of  the  exiting  jet.  The 

coefficient  B  was  empirically  determined  from  Jordinson  /Tl7 

2 

and  is  equal  to  0.19  (U^/O)  ,  where  Uj  is  the  jet  velocity 
and  U  is  the  mainstream  velocity.  To  further  describe  the 
problem  (see  Fig.  2),  a  system  of  natural  coordinates 
attached  to  the  jet  is  adopted.  Relative  to  the  fixed 
coordinate  system  and  an  arbitrary  point  [xp»ypJ  on  the 
surface  surrounding  the  jet  orifice,  the  natural  coordinates 
are  given  by 

K  ■  -  [zsina  +  (x-xp)  cosa J 

n  -  Yp  (5) 

C  ■  zcosa  -  (x-x  )  sina 

P 

where  €  is  the  coordinate  tangent  to  the  jet;  C,  the  coordi¬ 
nate  normal  to  the  jet  in  the  direction  of  the  center  of 
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curvature;  n*  the  coordinate  perpendicular  to  £  and  ;  ; 
and  a  ,  the  angle  between  the  x  and  £  directions. 

Another  extremely  important  part  of  Wooler's  formulation 
was  the  calculation  of  the  distribution  of  vorticity  within 
the  jet.  This  distribution  of  vorticity  was  determined  by 
adopting  methods  commonly  found  in  the  aerodynamic  theory  of 
lifting  surfaces  /—  227*  According  to  this  theory,  regions  of 
flow  external  to  the  jet  are  taken  to  be  irrotational  and  the 
deflection  of  the  jet  is  due  to  a  purely  inviscid  mechanism. 
This  inviscid  mechanism  is  expressed  as  a  balance  between  the 
forces  due  to  pressure  differences  across  an  element  of  the 
jet  and  the  centrifugal  forces  associated  with  the  jet  curva¬ 
ture.  These  pressure  differences  define  corresponding  velocity 
differences  according  to  Bernoulli's  Theorem.  These  resulting 
velocity  differences,  in  turn,  are  related  to  circulation 
according  to  Kelvin's  definition  of  circulation.  Therefore, 
with  the  above  assumptions,  the  circulation  is  related  to 
jet  curvature,  and  this  relationship  is 

tf-T-fO)2 

where  r  is  the  circulation  per  unit  length,  ds,  along  the 
jet;  R,  the  local  radius  of  curvature  of  the  jet;  and  d,  the 
diameter  of  the  jet  orifice.  Utilizing  the  empirical  tra¬ 
jectory  expression  /~Eq.  (4 )7  in  the  arc  length  derivative 
results  in 
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Incorporating  Eq.  (7)  into  Eq.  (6)  and  integrating  over  an 
element  of  jet  extending  from  to  Z2  results  in  the  follow¬ 
ing  expression  for  the  circulation  of  a  jet  element 

K  »  j  Ud  j^tan"1  e  2/Bd  -  tan-1  e  Zl/Bdj 

In  order  to  determine  the  flow  field  disturbance  on  the  plate , 
the  jet  was  divided  into  a  number  of  vortex  elements  whose 
strengths  were  calculated  using  Eq.  (8).  These  elements  are 
the  origins  of  characteristic  horseshoe  vortex  systems  where 
the  cross  member  is  bound  in  the  jet  and  the  trailing  sides 
are,  according  to  Wooler's  assumption,  tangent  to  the  jet  and 
separated  by  a  distance,  d,  the  diameter  of  the  jet  orifice 
(Fig.  3) .  Errors  associated  with  taking  the  trailing  vortices 
tangent  to  the  jet  instead  of  containing  them  within  the  jet 
as  in  the  actual  flow  situation  were  assumed  by  Wooler  to  have 
negligible  effect  on  the  velocity  field  along  the  plate.  Each 
of  these  horseshoe  vortex  systems  produces  an  interference 
velocity  on  the  surface  surrounding  the  jet  orifice.  To  de¬ 
termine  this  interference  velocity  at  any  arbitrary  point  on 
the  surface,  the  Law  of  Biot  and  Savart  £—2 37 


is  applied  to  each  bound  vortex  and  its  associated  trailing 
vortices.  For  the  bound  vortex,  -  d/2  <  L  ^  d/2,  and  for  the 
trailing  vortices,  <  L  <  ,  were  the  inteoration  limits. 

These  integrations  produced  closed-form  solutions  for  the 
interference  velocities  due  to  individual  jet  elements.  The 
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velocity  ratios  (Figs.  4/5,6)  show  fair  agreement  with  the 
experiments  of  Bradbury  and  Wood  /” 6_7.  However,  these 
plots  also  point  out  serious  deficiencies  in  the  model.  In 
the  far  field,  in  the  arc  35°  <  0  <  145°  (0  being  measured 
counter  clockwise  from  the  ray  extending  downs tream  from 
the  jet  origin) ,  the  method  represents  the  jet  interference 
effects  well,  but  in  both  the  upstream  and  downstream 
portions  of  the  surface,  there  is  a  definite  lack  of  agree¬ 
ment.  In  the  upstream  area,  the  blockage  effect  due  to  the 
presence  of  the  jet  in  the  uniform  flow  is  not  predicted, 
while  in  the  downstream  regions,  additional  wake  effects 
have  not  been  taken  into  account.  However,  in  view  of  the 
large  area  of  good  agreement  with  experiment,  Wooler's 
distributed  vorticity  model  was  selected  as  the  foundation 
for  further  investigation.  It  is  apparent  that  this  model 
gives  a  fair  representation  of  the  actual  contra-rotating 
vortices  associated  with  the  jet-crossflow  interaction. 
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IV.  EXTENSION  OF  WOOLER  MODEL 


Upon  successful  confirmation  and  correction  of  Wooler's 
distributed  vortex  model  and  as  a  step  towards  developing  an 
improved  model,  it  was  felt  necessary  to  verify  the  trailing 
vortex  assumption,  that  these  vortices  are  tangential  to  the 
jet  rather  than  contained  in  the  jet.  Retaining  the  basic 
theoretical  concepts  for  calculating  the  incremental  circula¬ 
tion  of  a  finite  length  of  the  jet,  containment  of  the  trailing 
vortices  was  approximated  by  successive  conforming  discretized 
vortex  elements.  These  vortex  elements  were  formed  by  restrict¬ 
ing  the  length  of  the  trailing  vortices  to  the  linear  distance 
between  z^  and  the  arbitrarily-chosen  end  points  of  the 
finite  length  of  jet,  ds.  In  order  to  satisfy  the  Helmholtz 
vortex  theorem,  that  a  vortex  cannot  end  in  the  fluid,  the 
above  horseshoe  vortex  system  was  closed  by  connecting  the 
trailing  arms  with  another  bound,  contra-rotating  vortex. 

(See  Fig.  7)  By  making  this  alteration,  vortex  loops  of  in¬ 
creasing  incremental  strength,  as  given  by  Eq.  (8),  can  be 
made  to  conform  to  the  curvature  of  the  jet,  as  in  the  actual 
flow  situation.  This  resulting  vortex  pattern  is  analogous 
to  that  created  by  a  continuously  accelerating  wing  ^"”23 _7. 

As  a  wing  accelerates,  the  strength  of  the  produced  vortices 
correspondingly  increase,  a  situation  very  similar  to  the 
change  in  elemental  vortex  strength  associated  with  the  jet 
curvature. 
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Using  this  contained-vortex  system,  the  interference 
velocity  due  to  the  jet  at  any  arbitrary  point  on  the  plate 
was  determined  by  applying  the  Biot-Savart  Law  (See 
Appendix  A).  In  this  application,  it  was  necessary  to 
approximate  the  trailing  vortices  as  straight-line  segments. 
However,  as  the  sizes  of  the  jet  elements  become  small,  the 
linear  approximation  can  be  made  to  conform  to  the  curva¬ 
ture  with  an  exactness  that  is  only  limited  by  numerical  prac¬ 
ticalities.  The  total  interference  components  due  to  the 
entire  jet  result  from  the  summation  of  elemental  contribu¬ 
tions  and  are  given  by 
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where  N  is  the  number  of  jet  elements  and  and  zJc+1  are 
the  endpoints  of  the  general  element.  Again,  the  parameters 
u,  v ,  and  w  are  defined  as  geometric  coefficients  resulting 
from  applying  the  Biot-Savart  Law  around  the  element.  The 
subscripts  1,2,3  and  4  refer  to  the  forward  and  after  bound 
and  left  and  right  trailing  vortex  contributions,  respectively. 
These  are  given  by 
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on  the  plate.  The  resulting  contours  of  constant  Dressure 
coefficient  are  shown  in  Pigs.  8,  9,  and  10. 

These  resulting  plots  show  a  somewhat  different  agree¬ 
ment  with  the  experimental  data  /*” 6  7  than  did  the  results 
of  Wooler's  model;  however,  the  discrepancy  can  be  directly 
attributed  to  the  manner  in  which  the  trailing  vortices  are 
treated.  In  Wooler's  model,  the  total  interference  velocity 
due  to  the  trailing  vortices  results  from  the  cumulative  el¬ 
emental  vortices  acting  at  different  distances  relative  to 
a  point  on  the  surface,  while  in  the  conforming  vortex  model 
the  trailing  vortex  contribution  results  from  the  cumulative 
elemental  vortex  strength  acting  at  a  single  radius,  (See 
Fig.  11).  Therefore,  Wooler's  assumption  of  treating  the 

trailing  vortices  as  being  tangent  to  the  trajectory  intro¬ 
duced  an  error  that  under  some  conditions  tended  to  improve 

agreement,  where  the  conforming  vortex  model  removed  this 
error  at  the  expense  of  some  loss  in  agreement,  at  least  at 
the  larger  values  of  U^/U.  However,  there  have  been 
studies  /“12,21 7  that  present  the  concept  of  an  effective 
origin  of  the  jet  vortex  system  -  a  region  where  the  turbu¬ 
lent  mixing  processes  extend  across  the  entire  jet  causing 
definite  deflection  of  the  jet  and  establishment  of  the  contra' 
rotating  vortices.  Schmidt  /”21 7  formulated  an  empirical 
expression  for  the  position  of  the  effective  source  as  a 
function  of  the  jet-to-main stream  velocity  ratio  given  by 

Za  -  (15  o/o)~l 
© 

where  o  *  U/U .  .  Applying  this  effective  source  expression 


to  the  conforming  vortex  model  produces  better  agreement 
with  experiment  (Fig.  12). 

Comparison  of  the  conforming  vortex  model,  including  an 
effective  source,  with  experimental  data  shows  the  same, 
if  not  better,  far-field  agreement  with  an  improvement  in 
the  predicted  pressure  distribution  in  upstream  and  downstream 
regions.  Therefore,  this  model  which  more  closely  represents 
the  actual  flow  situation,  replaced  Wooler's  model  as  the 
basis  for  further  investigation  and  improvement. 
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V.  IMPROVEMENT  OF  CONFORMING  VORTEX  MODEL 


Evaluation  of  the  conforming  vortex  model  results  in¬ 
dicated  the  continued  lack  of  agreement  with  experiment  in 
two  distinct  areas:  upstream  and  wake.  In  order  to  improve 
the  agreement  in  these  areas ,  it  was  felt  that  additional 
corrections  could  be  formulated  from  a  qualitative  analysis 
of  the  flow  phenomenon.  The  discrepancies  in  the  upstream 
area  can  be,  at  least  partially,  attributed  to  the  blockage 
of  the  mainstream  flow  by  the  presence  of  the  jet,  while  the 
discrepancies  in  the  wake  area  are  caused  by  turbulent  flow 
separation,  vortex  shedding  and  turbulent  entrainment.  In 
view  of  the  present  uncertainty  concerning  the  exact  character 
of  the  wake  area,  it  was  decided  that  the  development  of 
corrections  for  the  wake  must  await  future  investigation  and 
experimentation.  Therefore,  the  upstream  blockage  correction 
was  undertaken. 

Near  the  point  of  injection,  the  physical  appearance  of 
a  circular  jet  injected  into  a  crossing  flow  is  that  of  a 
cylinder.  Interpreting  the  jet  as  a  cylinder  in  an  irrotational 
flow,  potential  flow  theory  predicts  perturbations  to  the 
uniform  flow  that  are  given  by 


*»  ""  -y  COS  29 
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where  r  and  6  are  cylindrical  coordinates  of  points  in  the 
flow  field  and  a  is  the  radius  of  the  cylinder.  These 
velocity  perturbations  caused  by  blockage  can  be  combined 
with  the  interference  velocity  components  derived  from  the 
vorticity  formulation  to  adjust  the  coefficient  of  pressure 
(Eq.  9)  at  arbitrary  points  on  the  surface.  The  initial 
(and  most  simple)  assumption  was  to  make  the  diameter  of  the 
blockage  cylinder  the  same  as  the  diameter  of  the  jet  orifice 
The  resulting  pressure  distribution  indicated  an  excessive 
blockage  that  is  not  present  in  the  real  flow.  The  actual 
flow  does  not  produce  a  solid  interface,  as  would  a  cylinder, 
but  instead  an  entraining  interface  surrounding  a  core  of 
uniform  velocity  approximately  equal  to  U ^ ,  the  jet  exit 
velocity.  To  better  describe  this,  a  potential  cylinder  of 
diameter  less  than  the  diameter  of  the  jet  was  introduced. 
Using  many  combinations  of  cylinder  diameters  with  different 
jet-to-mainstream  velocity  ratios,  an  expression  for  a  in 
Eq.  (10)  was  formulated  in  terms  of  the  velocity  ratio. 

a  *  0.96  /a 

Applying  this  formulation  to  the  blockage-cylinder  perturba¬ 
tion  velocity  calculations  results  in  pressure  contour  plots 
(Fig.  13,14,15)  that  show  improved  agreement  with  experiment 
in  the  far  field  and  upstream  areas  for  jet-to-mainstream 
velocity  ratios,  4  <  i  <  8.  However,  as  the  velocity  ratio 
increases,  the  model  results  and  experiment  diverge.  The 
apparent  cause  of  this  disagreement  can  be  resolved  by 
qualitatively  examining  the  actual  flow.  With  increasing 


velocity  ratios,  the  entrainment  of  crossflow  fluid  by  the 
jet  becomes  the  more  dominant  interference  factor,  while  the 
blockage  factor  is  relatively  less  important.  Wu  and  Wright 
/T9_7  have  concluded  that  "the  amount  of  crosswind  fluid  en¬ 
trained  increases  with  increasing  jet  speed"  when  the  cross- 
wind  speed  is  held  constant.  In  addition,  this  increased 
entrainment  appears  to  be  especially  important  in  the  down- 
steam  wake  region  and  causes  a  further  reduction  in  the 
blockage  effect  as  compared  to  that  due  to  a  solid  cylinder 
in  a  crossflow.  The  resulting  effect  of  the  increased  rate 
of  entrainment  at  higher  jet  speeds  is  larger  interference 
velocity  perturbations  in  the  flow  field  surrounding  the  jet 
orifice  near  the  plate.  In  view  of  this,  the  lack  of  agree¬ 
ment  at  higher  velocity  ratios  can  be  explained  while  also 
nointing  out  the  need  for  additional  correction (s)  to  account 
for  the  change  in  the  relative  importance  of  blockage  and 
entrainment  with  changes  in  velocity  ratio. 
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VI.  CURRENT  MODEL  CAPABILITY 


Realizing  the  limitations  of  the  current  analytical 
model,  it  is  nevertheless  appropriate  to  investigate  its 
predictive  capability.  Figure  16  shows  the  results  of  the 
analysis  expressed  as  the  force  induced  on  a  flat  plate  due 
to  a  jet  injected  at  right  angles  to  a  crossflow.  The  force 
is  expressed  as  a  fraction  of  the  ideal  jet  thrust  and  the 
negative  values  indicate  suction  forces  (in  a  direction  op¬ 
posite  to  the  ideal  jet  thrust) .  Correlation  with  experiment 
/"" 4 _7  shows  a  difference  in  the  value  of  k  in  Eq.  (2).  The 
value  of  k  2  -4.4  resulting  from  this  analysis,  although  not 
in  the  predicted  range,  -3.0  <  k  £  -1.5,  still  displays  the 
correct  trend  and  provides  encouragement  as  to  the  potential 
success  of  a  model  based  upon  the  present  methods.  The  main 
difference  in  the  approximation  constant  results  from  the 
lack  of  proper  accounting  for  the  entire  entrainment  effect, 
a  factor  which  progressively  degrades  the  model  at  higher 
velocity  ratios.  However,  as  an  initial  design  approximation 
for  low  velocity  bow  thrusters  active  at  relatively  high 
ship  speeds,  the  conforming  vortex  model  presented  here  pro¬ 
vides  a  satisfactory  prediction  of  the  interference  effects 
—  a  step  toward  improving  the  accuracy  of  cowering  and  con¬ 
trol  determinations. 
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VII.  CONCLUSIONS 


The  conforming  vortex  model,  modified  with  corrections 
based  on  experimental  observations,  is  presented  as  an 
initial  formulation  of  a  simple  yet  physically  consistent 
representation  of  jet-crossflow  interference  effects.  Pres¬ 
sure  distributions  in  the  far- fie Id  are  adequately  predicted 
and  the  calculated  results  obtained  with  the  model  provide 
reasonably  good  agreement  with  experiment  for  the  lower 
jet-to-crossflow  velocity  ratios.  Agreement  is  less  than 
satisfactory  for  the  higher  velocity  ratios.  However,  the 
model  does  indicate  the  relative  importance  of  different 
flow  phenomena  with  changing  velocity  ratios.  The  most  evi¬ 
dent  weakness  of  the  model  is  the  lack  of  an  adequate  method 
for  including  entrainment  and  blockage  effects  near  to  the 
origin  of  the  jet  at  higher  velocity  ratios.  By  restricting 
the  range  of  application,  the  conforming  vortex  model  demon¬ 
strates  sufficient  agreement  to  be  utilized  as  a  first 
approximation  technique  for  designing  bow  thrusters. 


VIII.  RECOMMENDATIONS 


Further  experiments  should  be  conducted  to  provide  more 
correlation  data,  especially  in  the  case  of  water  jet  into 
water  crossflow.  Empirical  constants  utilized  in  the  current 
model  and  those  that  will  become  necessary  when  trying  to 
model  increased  entrainment  can  be  refined  when  a  larger  data 
base  is  available.  The  encouraging  results  of  this  method 
of  analyzing  the  flow  field  should  not  be  overlooked#  and 
additional  corrections#  as  indicated  in  this  analysis#  should 
be  made  to  improve  its  accuracy. 
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Figure  1.  Sketch  of  bow  thruster-crossflow  interaction 
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Description  of  coordinate  systems 


Comparison  Of  Pressure  Coefficient  Contours 
Experiment  - — /~6  7r  Conforming  Vortex  Model 
U./U  -  4.0.  “  “ 


Figure  9.  Comparison  Of  Pressure  Coefficient  Contours 
Experiment  —  /~6  7#  Conforming  Vortex  Model 
U,/U  -  8.0 
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Figure  10.  Comparison  Of  Pressure  Coefficient  Contours 
Experiment  — A  7»  Conforming  Vortex  Model 

Uj/U  ■  11.3 
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Figure  b.  Conforming  Vortex  Model 


Figure  11,  Comparison  Of  Wooler  (a)  And  Conforming  Vortex 

Model  (b)  Effective  Radii  For  Induced  Velocities 


Comparison  Of  Pressure  Coefficient  Contours 
Experiment  6_7r  Conforming  Vortex  Model 

With  Effective  Source  Correction — .U4/U  ■  8. 


Figure  13.  Comparison  Of  Pressure  Coefficient  Contours 

Experiment  —^”6  7,  Fully  Corrected  Conforming 
Vortex  Model — ,Uj/U  ■  4.0 


Figure  15.  Comparison  Of  Pressure  Coefficient  Contours 

Experiment  — /— 6  7»  Fully  Corrected  Conforming 
Vortex  Model — .0 .7u  ■  11.3 


APPENDIX  A 


FORMULATION  OF  INTERFERENCE  VELOCITY  COMPONENTS 

The  formulation  of  the  interference  velocity  due  to  the 
vortex  loop  system  of  the  conforming  vortex  model  is  given 
below.  The  formulation  requires  the  use  of  three  separate 
coordinate  systems  that  are  functions  of  the  endpoints  of  a 
jet  element  and  the  position  of  an  arbitrary  point  on  the 
plate.  From  Fig.  (Al) ,  these  are  given  by: 

5^  *  -  ^z^sino^  +  (x^  -  xp)  cosaj  (Al) 

n  -  yp 

»  Zjcoscij^  -  (x^  -  xp)  sinc^ 
for  the  upstream  bound  vortex 

53  “  ^z2sino2  +  (x2  -  Xp)  co3a2j  (A2) 

n  -  Yp 

C2  *  z2cosa2  -  (x2  -  Xp)  sina2 

for  the  downstream  bound  vortex ,  and 

5Bi  *  "  [Zi3inaB  +  (xi  "  Xp)  COSOIb] 

n  *  yp  (A3) 

?B  -  ziCosaB  -  (x.  -  xD)  sinaB 
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for  the  trailing  vortices,  where  x  and  y  are  the 

r  r 

coordinates  of  am  arbitrary  point  on  the  plate  and  x^,  z^ 

, z.;  are  the  coordinates  of  the  endpoints  of  an  element 
of  jet.  The  angles  a^,  a2  and  aQ  are  the  angles  between  the 
x-direction  and  respective  5's  and  are  given  by 

*  tan”1  [csch  (z^/B)j 

*  tan”1  [csch  (z2/B)j 

aB  *  tan”1  [(z2  -  z1)/(x2  -  x^] 

To  determine  the  induced  velocity  at  an  arbitrary  point, 
the  law  of  Biot  and  Savart 


L 


must  be  individually  applied  to  both  the  bound  and  trailing 
vorticies.  Due  to  similarities  in  the  integrations,  only 
one  bound  and  ine  trailing  vortex  formulation  will  be  pre¬ 
sented  and  necessary  changes  applied  for  the  other  vortices. 
The  circulation  of  a  finite  amount  of  jet,  given  by  /”Eq. 

(  8 )J 

w  /V\2  r  -i  z2/Bd  -i  vBdl 

K  *  j~  ud \jfLJ  [tan  1  e  i  -  tan  x  e  1  J 

is  the  incremental  contribution  which  is  added  to  each 
successive  element,  such  that  the  total  circulation  for 
a  given  element  is  given  by 
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Therefore#  the  induced  velocity  due  to  any  of  the  vortex 
filaments  associated  with  the  n-th  element#  1  £  n  z  N# 
where  N  is  the  number  of  jet  elements#  can  be  determined 
from 


(A6) 


where  L  is  the  length  of  the  filament#  8  is  the  angle 
between  the  vortex  filament  and  the  radius#  r,  to  an 
arbitrary  point  on  the  plate. 

For  the  forward  bound  vortex  {See  Fig.  (A2)) 
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Substituting  these  equations  into  Gq.  (A6)  gives 
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the  magnitude  of  the  induced  velocity  at  an  arbitrary  point 

on  the  plate  due  to  the  bound  vortex  of  the  n-th  element. 
The  x-direction  component  of  this  velocity  is  given  by 
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The  approximation  of  the  bound  vortices  as  being  linear 
eliminates  the  y-direction  component  of  induced  velocity 
and  since  an  image  system  i3  utilized  to  formulate  the  flat 
plate,  it  is  unnecessary  to  calculate  any  z-direction 
component  of  induced  velocity. 

For  the  left  trailing  vortex  filament,  as  viewed  from 
upstream  (See  Fig.  (A3) ) . 


Inserting  these  expressions  into  Eg.  (A6)  gives 


the  magnitude  of  the  induced  velocity  at  an  arbitrary 
point  on  the  plate  due  to  the  left  trailing  vortex  of  the 
n-th  element.  Resolving  this  velocity  into  x  and  y  com¬ 
ponents  gives 
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Inserting  Eq.  (A12)  into  Eq.  (A13)  and  applying  integration 
limits  gives 


Integration  of  Eq.  (A14)  yields 


(A15) 


the  x  and  y  comoonents  of  induced  velocity  due  to  the  left 
trailing  vortex.  Similarly  the  x  and  y  components  of  the 
induced  velocity  due  to  the  right  trailing  vortex  are  given 
by 
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(A16) 


Therefore,  the  total  x  and  y  components  of  induced  velocity 
due  to  an  entire  element  of  jet  are 

qxT  *  qxl  +  qx2  +  qx3  +  qx4 

(A17) 

qyT  "  qy3  +  qy4 

Inserting  Eqs. (A10) ,  (All),  (A15),  and  (A16)  into  Eq.  (A14) 
yields 

Fn  r 

qxT  *  TiF  [U1  cosai  "  u2  coso2-w1sina1+w2sino2+  (-w3+w4 )  sinaB] 
r  (A18) 

V  *  -o  [v3  -  v<] 
where 
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“i  -  a  **  i\  n  * 
«i  *  ^  LvTT— t; 


kl^lW*)2  vS12+?12-*'(n-si) 


n  -  *  ] 

u«.  2 .  i _ iv  2  J 


J2 _ J  n  + 


?22+C22+(n+%)2  Vc22+c22+(n-»j) 


i—l 

2  J 


)  -  h  r _ sb2 


j  ■«I1+'b2+  "’-',i  2 


1  +  S  J  _ 

♦  "’+,»2Lv£L«.2+ 


*Bl  1 


CB  +  (H+%)  Ly522+i;B2+(n+H)2  V?Bl+^B2+(n+ls)2j 

Substituting  the  expression  for  I*  from  Eq. (A5)  into  Eqs. 
(Al8 )  and  non-dimensionalizing  yields 


^(r)  Y  [“wl8inVw28ina: 


i+l+  f”w3+w4)  sina^+UjCOs  a^- 


UjCoaai+1  x  ^  tan”1  f  * 

^  i _ i  \  i 


8k+l  *k 
i  -  e 


1  +  e 


‘k+l+’k 


v  _  1 
U  S 


/  tlA  2  N  f  i 

(<r)  t '  [v3  -  v<]  *  t 


1  +  e 


Vrt! 
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THIS  PROGRAM  SIMULATES  A 
FLUID  JET  INJECTEC  INTO 
A  UNIFORM  CROSS  FLO  I*  ANO 
FURTHER  DETERMINES  THE 
PRESSURE  01 STR I  ALT  I C  N  ON 
THE  SURFACE  FROM  WHICH 
THE  JET  IS  EJECTEC 


X 

X 

X 

X 

X 

X 

X 

X 

X 


xxxxxxxxxxxxxxxxxxxxxxxxxxxx 


DEFINITION  CF  TERMS 


NTYFE  -  MODEL  TYPE JC-CCNFQRM 

1- CCNFORM  BLOCKAGE 

2- RAISEO  CONFORM 

3- RAISED  CONFORM  BLOCKAGE 

(INPUT) 


C 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NR  AC 

NANG 

NPCINT 

LJ 

LM 

C 

P 

THETA 

>P 

YP 

X 
2 

21 NC 
CELTA2 
CELCCK 
HKSJV 

XI 

iET  A 
ETA 

ALPHA 

FAREA 


NUMBER  OF  RACIAL  PQSITICNS  CN  SURFACE 
( INPUT) 

NUM8ER  OF  ANGULAR  POSITIONS  ON  SURFACE 
( INPUT) 

NUMBER  OF  ARBITRARY  POINTS  CN  SURFACE 
( NRAO+NANG ) 

JET  SELOCITY  UNPUT) 

MAIN  STREAM  VELOCITY  (INPUT) 

ACTUAL  DIAMETER  OF  JET  CRIFICE  (INPUT) 

NON-CIMENS ICNAL  RAOIAL  POSITION  OF  SURFACE 
POINTS  (DIMENSIONED  AT  LEAST  NRAO)  (INPUT) 

ANGULAR  POSITION  OF  SURFACE  POINTS 
(DIMENSIONED  AT  LEAST  NANG) 

X-POSITION  CF  SURFACE  POINTS 
(DIMENSIONED  AT  LEAST  NPCINT) 

Y-POSITICN  CF  SURFACE  PCIMS 
(DIMENSIONED  AT  LEAST  NPCINT) 

HCRI20NTAL  CGCRDINATE  OF  JET 

VERTICAL  COORDINATE  OF  JET 

NON-CIMENS ICNAL  INCREMENTAL  STEP  SIZE  FCR  Z 

2-PCSITION  CF  EFFECTIVE/VIRTUAL  SOURCE 

EFFECTIVE  BLOCKAGE  CYLINDER  Cl  AM  ET  E  P 

CUMULATIVE  VCRTICITY  ALCNG  THE  JET 

NATURAL  CCCPCINATE  TANGENT  TO  JET 

NATURAL  CCCRCINATE  NORMAL  TC  JET 


NATURAL  COORDINATE  PERPENC 1CULAR  TC 
ANO  2ETA 

ANGLE  BETWEEN  X-AXIS  ANO  XI-AXIS 
INCREMENTAL  PLATE  AREA 


>1 
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TCTFCR  -  TOTAL  RESULTANT  FORCE  ON  PLATE 
>eAP  -  X-PGSITION  CF  RESULTANT  FCFCE 
VT  -  NORN ALIZEC  INCUCED  FORCE 
IMPLICIT  REAL<8  (A-P,C-Z> 

^CI PENSION  UISC^SOJ.VISO.SO.WISOtSCJ.RISO)  ,ThETA(50>, 

*CIPENSICN  Xp|«5&)?YP  (450  *ZP( 450  i 

DATA  INPUT 

REAC(2t60  )  NR  AD 
PEAC  (2*601  NANG 

ZERO  CUT  STORAGE  MATRICES 

NPCINT  «  NR AQYNANG 

CC  1C  I  »  I fNFCINT 
XP(IJ  »  0.0 
YP< I)  -  0.0 
ZP(I)  *  0.0 
CONTINUE 

CC  30  I  -  I »NRAO 
RCII  *  0.0 
CO  2C  J  »  I, NANG 
THETA(J)  «  O.C 
UltJI  -  0.0 
VII.J)  *  C.C 
CP(I,J)  «  0.0 
CONTINUE 
CONTINUE 

REAC(2»60)  NTYPE 

REACJ2.40I  UJ,UP,0 
NR  ITE  (6.50)  U«{ »UMf  D 

REAC(2»70I  ZINC 

RE  AC (2t70 ) (R( I)«I>lfNRAC) 


VELOCITY  RATIO 
VR  «  LJ/UM 

PC SIT I CN  OF  EFFECTIVE  SOURCE 
CELTAZ  ■  VR*0SQRT (VRJ/L5.0 
CORRECTION  FOR  0ELTA2  >  0.0 
IF<NTYPE.LE.II  CELTAZ  ■  0.001 


NRITE(6'8CJ  OELTAZ 
NR  I TE  (6  1 90)  ZINC 

£Q5^T(3F8.4) 

5C  „F0PPAT('0»t2X,»JST  VELOCITY  -  »,  F8 .4,2X ♦ • UN  I FCRP* , 
X  'VELOCITY  ■  •  ,F8.4,2X»  'JET  Cl AMETER  *  «,F8.4» 


FCRP  AT (141 
FCRNAT(F1 C. 5) 

FCPPAT(«0,,2X,'EFFECTIVE  SO 


FCRPAT(«0',2X1«EFFECTIVE  SOURCE  IS»,F 


IS*  *F5.2v  2X» 


DEFINITION  CF  PI 
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5 

C 


PI  »  2.0*ARSINU.00CCC000) 
TRAJECTORY  CONSTANT 
E  »  C.1S*VR**2 

CALCULATION  OF  INTERFERENCE  VELOCITY 
L  »  C 

CALCON  ■  OFLOAT(NANG)  -  1.0 

INCREMENT  RAO  IU  S  (INPLT  VALUES) 

CC  140  I  *  l'NRAD 
RCALC  =  R ( I ) 

INCREMENT  THETA  (5  OEC-REE  INCREMENTS) 
1*  NANG 

DFLOAT(K)*  ( PI/CALCON  ) 
SHIFT  TO  CARTESIAN  CCCRCINATES 


CO  130  J 
K  *  J-i 
THETA(J) 


l  »  L+l 

COST  »  OCOS(THETA( J ) ) 
SINT  »  OS  IN (THETA! J ) ) 
XP  (  L  I  »  RCALC  *COST 
YP 4  L  )  *  RCALC*SINT 


11 

22 


0.0 

C.00CG1 


ICC 


HKSAV  *  0.0 
Z2 


21 

2 


11  *  Z] 
22  «  12 


Z 1  ♦  ZINC 


Oi 

Oi 


LTAZ 

LTAZ 


>1 

>2 


B*(DC0SH(Z1/B)-1.0) 

Ij-i. 


«  B*(0C0SH(Z2/Bi 


0) 


TANAB  «  ( Z2-Z1)/(X2-X1) 

TANAd  *  1.0/0SINHCZ1/B) 

TANA2  *  l  .0/0SINH4  Z2/E) 

ALFCNE  *  CATAN(TANAl) 

ALFTHO  «  CATAMTANA2 ) 

ALFEAR  «  OATAN (TANAS ) 

SINAI  •  DSIN(ALFONE) 

| IN A2  «  DSI N ( ALFTWG ) 

SINAE  *  CSIN(ALFBAR) 

COS/1  «  CCCS(ALFONE) 

15^2  «  OCOS ( ALFTWG ) 

lAe  »  occs(alfbar) 

IN  TERMS  OF  CARTESIAN  COORDINATES 


NATURAL  COORDINATES 
X 


liftt 


il 


BARI 

E/R2 


-(Z1HSINA1 
-I Z22*SINA2 


♦ 

-(Zll*SINAE  ♦ 
-(222*SINAE  ♦ 


(X1-XP(L)I*CCSA1) 
<X2-XP(L))*CC$A2) 
(Xl-XP(L) >*CCSABJ 
(X2-XP  (L)HCCSAB) 


ETA  «  YP4L) 

2ETA1  ■  (Zll*COSAi  -  ( Xl-XPCL ) ) *SIKA 1 ) 

if})? :  liff.*{?!ii :  lirtfitiCiltlli 
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C  GEOMETRIC  CONSTAMS 


ETPl 
ETMI 
SA  * 
SE  * 

Sc  « 

SC  * 
SE  « 
«F  * 

Sg  » 

SI*  « 
SI  « 
C  A  * 
C8  * 

cc  * 

CE  * 
CF  * 
CG  * 
CH  * 
Cl  « 

CCM 

CCN2 

COM3 

CQN4 


»  ETA  «  C.5 
»  ETA  -  C.5 
XICNE**2 
XITW0**2 
X IBARl**2 
XI 8AR2**2 
2ETA 1**2 
ZETA2**2 
2ETAE**2 
ETPL**2 
ETMI**2 
OSCRT ( SA  ♦  SE 
DSQRTl SA  ♦  SE 
OSCRTISe  ♦  SF 
DSQRT(  SB  ♦  SF 
OSQRTl SO  +  SG 
CSCRT1SC  ♦  SG 
DSQRT (SC  +  SG 
CSCRT( SC  ♦  SG 


SH) 

SI ) 
SH  ) 
SI  ) 
SI ) 
SI  ) 
SH  ) 
SH) 


(ETPL/CA  -  E7HI/CB)/(SA  ♦  SE) 
(ETPL/CC  -  E7PI/DE)/(Se  ♦  SF) 
(XIBAR2/OF  -  XI EARi/CG)/ (SG  *  SI) 
( X I8AR2/0H  -  XIBAR1/OI ) / ( SG  *  SH) 


LI  «  ZETA  1*CCM 
12  *  ZETA2*CCN2 
Vi  *  ZETAE*CCN3 
V*  «  ZETA0*CON4 
H  *  XICN£*CGM 
V> 2  «  XITWC*CON2 
V3  *  ETMI *CQN  * 
V4  ■  ETFL*CCNA 


C 

C 


C 


C 

8 

C 


INCREMENTAL  VQRTICITY  TERM 

FF  »  ((DEXP(Z2/B)-0EXP(Z1/B))/(I.C-»CEXP((Z2  +  Z1)/E))) 
FK  »  41. 0/8*0)*  VR**2  *  DAT  AN (FF  ) 

FK  «  HK  ♦  HKSAV 
FKSAV  »  HK 

INCREMENTAL  INTERFERENCE  VELOCITY  COMPONENTS 
LP  «  C.O 

LP  »  FK*((-mASINAl+N2*SINA2+(-W3+M4)*SINAB+Ll*CCSAl) 
X  **  ( L2*COSA2 ) ) 

VP  *  0.0 

VP  «  HK*(V3  -  V4) 


U(l'J)  *  UP 
»J)  ♦  VP 


V(  I 

TEST  FCR  CONVERGENCE 


TEST  «  CABS(l)F)  *  DABS(VP) 


IF  (TEST . LE.C.COOOCI)  GC  TC  liO 
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■a 


c 

110 


GC  1C  100 


CCNTIMJE 

CALCUL AT ION  OP  BLXKAGE  VELOCITY  TERN 
IFINTYPE.EQ.C )  GO  TC  120 
IFINTYPE.EC.2  )  GO  TC  120 
PeiCCK  '■  C.96/CSQRT  IVR  ) 

C8LCCK  *  2.0*PBL0CK*C 


PK  »  (IRBLOCK/RI I) )**2)/2.0 
LC  *  RK*(SINT<*2  -  CCST**2) 
VC  *  -2.0*ftK*SINT*CCST 

U  I,JI  *  U(  I, J)  *  uc 

V(I f J)  ■  VI  I.J)  ♦  VC 


C  CALCLL AT  ION  CF  PPSSSURE  CCEFFICI ENT 
120  CPIItJI  »  -4.C*(UCI,J)  ♦  U  ( I  tJ  1**2  *  V(I,J)M2) 

C  CFANGE  CES IGNAT ICN  OF  VARIAELE 

2 P  <  L  I  »  C  F<  1 1  J  ) 


C 

120 

140 


CONTINUE 

CONTINUE 


C  CALCULATION  OF  INCREMENTAL  FORCE 


RG  «  0.0 
PI  *  0.0 
AREA  »  0.0 
PAREA  *  0.0 
FCPCE  »  0.0 
YMGN  ■  0.0 
TCTFCR  *  0.0 
TYNCN  *  0.0 


NAREA 


NFAC  -  1 


CC  170  I«1,NAFEA 
K  *  I  ♦  1 
PI  *  PCI) 

RG  *  R(K ) 


C  CALCULATION  CF  INCREMENTAL  AREA 


AREA  ■  PI*( R0**2  -  P  1**2) 
FAREA  *  AREA/ <2.0*CALC0NI 
CO  160  J*l» NANG 


C  CALCUIATICN  CF  AVERAGE  PRESSURE  COEFFICIENT 


CPASG  «  (CPdtJl  ♦  CPI  K  t  J  H  /2  .0 
CALCULATION  Of  EFFECTIVE  FCRCE 
FORCE  «  PAREAACPAVG 


IFI..EQ.1)  GC  TO  15C 
IF  (J.EQ.NANG)  GO  TC  ISO 
FCRCE  «  2.0AFCRCE 


C  CALCUIATICN  CF  MCMENT  ARM 
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15C  X^  *  ((PC  ♦  RI)/2.C)*CCCS(THETA(JI  ) 
C  CALCULATION  CF  MOMENT  ABOUT  Y-AXIS 
C ►  »  FOPCE*XM 

C  SUMMATION  CF  TOTAL  FORCE 

7CTFCR  *  TOTFCR  ♦  FCRCE 
C  SUMMATION  CF  TOTAL  MOMENT 
TYMCM  a  TYMCM  ♦  YMCM 


160 

UC 


CONTINUE 
CCMINLE 

CALCLL ATICN  CF  CENTER  CF  ACTION 

>E  AF  »  TYMCM/TOTFQR 

CALCLLATICN  CF  NORMALIZED  INDUCED  FORCES  Y/T 

\RINV  «  1.0/VP 

YT  *  2.0*TCTFCR*VRINV**2/PI 

DATA  CLTPUT 

WRITE  (6.260 )  C0LOCK 
WRITE<6,220)  TCTFOR.XEAR 
WRI TE (6. 2 AO  I  YT»  VRINV 


CC  220  I  »  l  •  N  AN  G 
angle.*  TtETM.U*iac.o/Pl 

«,F6.2.8X.'R'.9X.'CP».C.X.'L'» 


180 

WRITE(6,190i  ANGLE 
19C  FGRMATCO' »2X  .'THETA 
X  9X »  *  V  ) 

CC  <10  J  »  1 .  NR  AO 

WRITE (6, 2 COl  R(J).CF(J.  1 1  *  U ( J  » I)  * V ( J  •  I ) 

2CC  FORMAT! '  '  ,23X,F6.2»2X,F8.2,3X,F6.2.AX,F6.2) 

210  CCNTINUE 
220  CONTINUE 

230  FCRMAT('0»,2X, 'RESULTANT  FORCE  OF  •  , IX, F7. 2 . 1 X  , 

X  •  ACTS  AT  X  a  »  ,F6.2  I 

2A0  FGRMAT('0».2X,'Y/T  a  »  ,F6 .2 ,1 X»  'FCR  LM/UJ  »  •  ,F6.3> 
WRITE  (7*2  5CMXP(  II  tYPUI.ZPU  I  . 1*1 ,  NFOINT I 
250  FORMAT! 3F20  .A  ) 

260  FCFWAT('0«f2X. 'EFFECTIVE  BLOCKAGE  CYLINDER* . 

X'CIAfETER  *  '  .F4.2) 

27C  STOP 
ENC 
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